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Abstract 
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1 Introduction 



The nonequilibrium dynamics of quantum systems are a timely and relevant 
subject that is important to many areas of study in physical science, from con- 
densed matter systems to cosmology. In most cases we are usually interested 
in the processes of thermalization and equilibration, the determination of the 
relevant time scales involved, together with an understanding of the genera- 
tion of entropy and particle production in nonequilibrium dissipative systems 
interacting with an environment. However, despite the importance associate 
to these processes, nonequilibrium problems are still poorly understood (see, 
e.g. [1]). The nontrivial nonequilibrium dynamics of fields, for instance, have 
diverse applications, finding use e.g. in the studies concerning the recent exper- 
iments in ultra-relativistic heavy- ion collision [2]; applications to the current 
problems of parametric resonance and particle production in cosmology [3] ; or 
in the context of the recent studies involving the intrinsic dissipative nature 
of interacting fields [4,5,6]. In addition to that, typical problems we have in 
mind to study are those related to the nontrivial out-of-thermal equilibrium 
dynamics associated with phase transitions in different physical systems. As 
a few examples we may cite include the current applications to the study of 
formation of Bose-Einstein condensates after a temperature quench [7] , or in 
the study of the dynamics of coupled fields displaced from their ground states 
as determined by their free energy densities [8]. 

One problem that we will be particularly interested in here is that one that 
involves quenched phase transitions, either caused by a sudden change in tem- 
perature of the system, or by a time change of the parameters of the model 
(like frequency or couplings), in which case we have to deal with situations 
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characterized by mixed states involving general density matrices and where an 
equilibrium or quasi-equilibrium formalism do not apply. Common tools used 
to study such nonequilibrium quantum systems are the Schwinger's closed 
time path functional formalism and the canonical quantization methods [9]. 
In this paper we will mostly be concerned with the understanding of simple, 
exactly soluble quench scenarios from the stand point of the Liouville-von 
Neumann (LvN) quantum canonical formalism [10]. In this formalism we have 
that the time evolution of the density operator p is given by the LvN equation 
with Hamiltonian operator H: 

^p(t)=i[p(t),H] (1) 

and satisfying the usual normalization condition Trp = 1. In order to solve 
Eq. (1) we can make use of the following fact. Suppose any two operators 6i(t) 
and o 2 (t) that satisfy the LvN equation (j = 1,2) 

^d j (t) = i[d j (t),A\, (2) 

then we can easily show that any operator function of 6i(t) and 62 (i), /(oi, 62), 
will also satisfy Eq. (2) provided that the time dependence in / only comes 
through its dependence in 61 and 6 2 . This important property allows us to con- 
struct an exact solution to Eq. (1) as a function of time dependent operators 
that satisfy the LvN equations. The LvN approach has been recently used in 
Refs. [11,12] to treat different problems in quantum field theory and cosmol- 
ogy and also emphasizing the advantages of the method compared to other 
approaches for nonequilibrium dynamics of quantum systems. The method 
may also be more efficient than the more common functional integral method, 
where most of the time we have to deal with nonlocal terms in time that can 
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only be appropriated deal with by recurring to linear response and/or adia- 
batic approximations (by requiring that the system moves sufficiently slowly 
as compared to typical microscopic time scales) that may not always apply to 
physically interesting systems (see for example Refs. [4,5,6]). 

In this paper we will study the use and applicability of the LvN approach for 
solving different systems of quantum oscillators subject to different quenches 
that are exactly soluble. With the study of these systems we expect to gain 
enough insight to then extend this analysis to more complex problems, in- 
cluding the nontrivial cases involving nonlinear interactions, where usually we 
must make use of perturbative methods, like in quantum fields in interaction 
and in the description of system evolution (like through a time dependent or- 
der parameter), when interaction with a thermal environment is included. In 
the next section we will give a simple illustration of use of the LvN approach 
by solving the (frequency) time dependent harmonic oscillator. Then in Sec. 
Ill we will study two examples of interacting quantum systems composed of 
two harmonic oscillators, initially in different thermal states, that are then 
coupled linearly at some given time. The case of a smooth quench is explic- 
itly solved. Then in Sec. IV we study the case of a system, represent by an 
harmonic oscillator, coupled linearly to a bath, represented by an infinity set 
of other oscillators, like in Caldeira-Leggett models and we explicitly evaluate 
the number density for the system, in the case of a sudden quench, in the 
context of the LvN approach. In Sec. V we discuss a general result concerning 
the addition of nonlinear interactions, for a model with quartic interactions, in 
the context of the approach discussed in the previous sections. Finally in Sec. 
VI we give our concluding remarks. Two Appendixes are included to clarify a 
few technical details. 
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2 The Time Dependent Harmonic Oscillator 



Let us begin by showing how the LvN approach works in the simplest prob- 
lem of the time dependent harmonic oscillator case, described by the time 
dependent Hamiltonian: 

h = \f + (3) 

with a time dependent frequency u{t). We consider that initially the system 
is in thermal equilibrium state at a temperature T = 1/(5 (note that for 
convenience we use throughout this paper the natural units, ks — h — 1). 
Next, we introduce time dependent annihilation and creation operators, a(t) 
and at(t), given by 

a(t) = A(t)q + B(t)p (4) 

and 

a\t) = A*(t)q + B*(t)p, (5) 

where the time dependence comes through the c- number coefficients A(t) and 
B{t). From the relation 

[a(t),a\t)] = l, (6) 
we obtain that 

A(t)B*(t) - B(t)A*(t) = -i. (7) 
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Using Eq. (4) in Eq. (2), we obtain 

A(t)q + B(t)p = -A(t)p + cu 2 (t)B(t)q, 



(8) 



where dots mean time derivatives. From Eq. (8) we then get that 

B{t) + u 2 {t)B{t) = (9) 

and A(t) is related to B(t) by the relation 

A(t) = -Bit). (10) 

Using the above equation we can rewrite Eqs. (4), (5) and (7) and express 
them as 

a(t) = B(t)p-B(t)q, (11) 

a f (f) = B*(t)p- B*(t)q (12) 

and 

B(t)B*(t) - B(t)B*(t) =i. (13) 

It is easy to see that Eq. (13) is the Wronskian of Eq. (9), i.e., warranting 
that the operators a(t) and ^(t) are annihilation and creations operators at 
all times. Now since our annihilation and creation operators satisfy the LvN 
equation, any function of them must also equally satisfy it too. In particular 
we can construct an exact solution for p(t). Since our initial condition is one 
of thermal equilibrium we can, therefore, write 

p(t) = ^- exp j-/foo ^(t)a(t) + 1] J , (14) 
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where 



Zp = Tr exp |-/3w 0^)6(0 + ^ } 



(15) 



and uo = cu(t — > — oo), z.e, we are choosing as initial condition that our system 
is at thermal equilibrium at t — > — oo with temperature 1/(3. Since relation (6) 
is satisfied, we can define a time dependent number operator N(t) = a^(t)a(t), 
with time dependent number eigenvectors \n(t)). However, because Eq. (6) is 
a time independent relation, the number eigenvalues n are time independent, 
i.e, we have 



N(t)\n(t)) = n\n(t)). 



(16) 



Using Eq. (16) it is easy to show that Zp is the partition function for the 
initial thermal equilibrium state, 



Z P = J2 ( n W I ex P \ -fao 
= gexp{-^ (n + i)} 



\n(t)) 



2 sinh ( 



Pup 
2 



(17) 



Note that p(t), given by Eq. (14), solves exactly the out-of-equilibrium dynam- 
ics of the system with the condition of thermal equilibrium as initial state. In 
order to compute the time dependent thermal averages of certain operators it 
is convenient to rewrite Eqs. (11) and (12) in terms of the initial annihilation 
and creation operators, that we will denote simple as a and a'. For this end 
we first solve for q and p in terms of the initial values of B(t) and B(t) that 
we will denote just as B and B, respectively, without the time argument. We 
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then find that 



q = i(B*a- Bat), 
p = i(B*a- Ba r ). 



(18) 
(19) 



Using now Eqs. (18) and (19) in Eqs. (11) and (12), we then obtain the Bo- 
goliubov transformation: 



a(t) 



aUt) 



I \ 

a(t) /3(f) 



P*{t) a*(t) 



( \ 

a 



\ J 



(20) 



where a(t) and P(t), the Bogoliubov coefficients, are given by 



a(t) = i (B*B(t) - B*B(tj) , 
(3(t) = i (BB(t) - BB(t)) . 



(21) 
(22) 



In terms of ct(t) and f3(i) Eq. (13) gives then the following relation: 



l«(t)| 2 -|/3(t)| 2 = l. 



(23) 



The utility of Eq. (20) is clear once we compute the thermal expectation value 
of the number operator N, 



n(f) = Tr[p(f> 

l 00 r 

Z /3 „=0 1 



a ] {t)a{t) + 



a) a \n(t)). 



(24) 



Solving Eq. (20) for a and a* in terms of a(t) and at(t), substituting them in 
Eq. (24) and using Eq. (16), it is easy to find that 

n(t)= \\a(t)\ 2 + \P(t)\ 2 }n+\P(t)\ 2 , (25) 



where 



n = -=-i (26) 



is the initial thermal equilibrium average for the number operator. 
2. 1 The Time Independent Harmonic Oscillator 

Let us first consider the trivial case of the time independent oscillator, with 
ui 2 (t) = uj 2 as an example of application. In this case the solution for Eq. (9) 
is just 

B{t) = ae iu}ot + be-™ ot . (27) 

substituting the above expression for B(t) in Eq. (13) we obtain 

= (28) 

In this case we can obtain a time independent number operator choosing b = 
and \a\ = 1/^2uj . The other choice a = and |6| ^ is simply ruled out 
since in this case we would get \b\ 2 < 0. Choosing a null phase for a, we then 
obtain the solutions 

a(t) = -j={p-iu Q q) (29) 
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and 



-iiOQt 



a\t) = -j== (p + iuj q) . 



(30) 



Although the creation and annihilation operators are time dependent the num- 
ber operator N = at(t)a(t) is time independent, as expected in this trivial 
example. 

2.2 An Exactly Soluble Time Dependent Harmonic Oscillator 



Let us consider now the case of an harmonic oscillator with the time dependent 
frequency: 

,2 



cosh (t/r] 



(31) 



Here the function B(t) satisfies the differential equation 



B(t) + 



cosh 2 (t/r) 



B(t) = 0. 



(32) 



Eq. (32) can be solved with the condition that at t — > — oo the system is 
at thermal equilibrium at temperature T. From the results of the previous 
subsection we then have to solve Eq. (32) with the asymptotic condition that 

iu t 



B{t -> -oo) 



(33) 



The above problem is formally identical as an one-dimensional quantum scat- 
tering problem. The solution of Eq. (32), under the asymptotic condition given 
by Eq. (33) is [13] 



B(t) = -j= [sech(t/r)r° F 



a,b, c, i + -tanh(t/r) 



(34) 
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where F(x, y, z, w) is the Gauss' hypergeometric function [14] and the coeffi- 
cients a, b and c are given by 



1 1 

a = - + iujqt + Ju 2 t 2 + -, (35) 



fe = i + ^ r-X 2 r 2 + i (36) 



and 

c—l + iu r. (37) 

When t — > oo intuitively one expects that the system will thermalize. Taking 
the t — > oo limit in Eq. (34) we obtain the result 

m >oq) r(~c)T(~a + b-c) e*^ T(c)T(c — a — b) e~ iu}Qt 
{ ' T(a)T(b) V^~o T(5-~a)T(5-b)VM>' 

Substituting Eq. (38) in Eqs. (21) and (22) we obtain for n(t), given by Eq. 
(25), an expression with oscillatory behavior in general. For these cases clearly 
the system never reaches thermal equilibrium. But if we analyze with some 
care the arguments of the gamma- functions in Eq. (38), we see that for those 
values of r such that r 2 u 2 = N(N + 1), with iV integer, then the second term 
in Eq. (38) vanishes. In this case the system reaches a final thermal equilib- 
rium state with final temperature that can be shown to be the same as of the 
initial one. This same result was also obtained previously in Ref. [15] whose 
authors have studied this problem solving the LvN equation obtained from 
a variational principle. We can also show that in the final thermal equilib- 
rium state there will be no dependences on parameters other than those of 
the initial thermal equilibrium state (c^o and (5). All the factors that appear 
multiplying the exponential factor in the first term in Eq. (38) contribute 



11 



with a time independent phase factor. This is equivalent to what happens in 
the one dimensional scattering problem as consequence of the conservation of 
probability in a reflectionless potential [13]. In computing physical quantities 
the phase factor cancels out and then, when the system reaches a final thermal 
equilibrium state it loses all the information about how this state was reached. 
Implicitly we have here a manifestation of decoherence already at the level of 
this simple problem as a consequence of a particular quench scenario and that 
will be so important when models with larger degrees of freedom are taken 
into account. 

As an example of a concrete result, let us consider for convenience the case 
N — 1, i.e, when t 2 uj 2 = 2. In this case B(t) reduces to 

pivot 
y/2uj (l + IUJqT) 



Substituting Eq. (39) in Eqs. (21) and (22) we get the results 



Ht)p 



1 - tanh 2 (t/r) 



4u 2 r 2 (l + u 2 T 2 ) 



(40) 



and 



1 - tanh 2 (t/r) 



4uj 2 t 2 (1+u 2 t 2 ) 



(41) 



Using these equations for |a(t)| 2 and |/3(t)| 2 in Eq. (25) we then obtain for 
n(t) the result: 



n(t) =11 + 



1 -tanh 2 (t/r) 



2u 2 t 2 (1+u 2 t 2 ) 



n + 



1 - tanh 2 (t/r) 



4u 2 t 2 (1+u 2 t 2 ) ' 



(42) 
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Results for n(t) for u (3 = 1 and for different values of the "decay time" r are 
shown in Fig. 1. Note that the peak at t = just represents the point around 
which the time-dependent frequency u(t), Eq. (31), reaches its maximum value 
and the peak is higher the smaller is the decay time, in agreement to the fact 
that, from Eq. (41), the smallest is r the largest is the production of excited 
states (particles, in the context of quantum fields), as represented by |/3(t)| 2 . 



3 Time Dependent Quantum Interacting Oscillators 



Lets us now consider a system of two time dependent oscillators, linearly 
coupled with interaction strength X(t), with Hamiltonian given by 

2 

2 



(43) 



J'=l 



The time dependent annihilation and creation operators can now be written 
as (i = 1,2) 

2 



(44) 



and 



j'=i 



(45) 



In an identical way as we have performed in the single time dependent oscil- 
lator of last section, we require that the operators a,i(t) and a\(t) satisfy the 
LvN equation, Eq. (2). From this requirement the following equations for the 
time dependent coefficients B^t) follow: 



B^t) + uKQBait) + X(t)B i2 (t) = 



(46) 
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and 



B i2 (t) + u%(t)B i2 (t) + \(t)B n (t) = 0. (47) 
The condition 

[a l (t),a](t)]=5 lj , (48) 

leads now to the equation 

£ [B ik (t)B* jk (t) - B ik (t)B* jk (t)} = ib l3 (49) 
k=i 

and the identities 

[ai(t),aj(t)] = [a\(t),a](t)]=0 (50) 
lead, respectively, to the following additional equations for B^t): 

2 

J2 [B ik {t)B jk (t) - B ik (t)B jk (t)] =0, 
k=i 

E [B* k (t)B* k (t) - Bt k (t)B] k (t)\ = 0. (51) 
k=i 

Using Eqs. (46) and (47), it can be shown that Eqs. (49), (50) and (51) must 
be time independent. As in the one single harmonic oscillator case we can 
write the time dependent annihilation and creation operators as functions of 
the initial annihilation and creation operators (that we just call by and a\, 
respectively, without the time argument, as indicating that these quantities 
should be evaluated at the initial time). This can be accomplished by using 
the requirement that initially the system consists of two independent harmonic 
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oscillators with frequencies u\ and 0J2, each one at well defined thermal equi- 
librium states with temperatures and I/P2, respectively. In this case the 
matrix {B}, at the initial time, has elements B^ given by 



B{j Sij - 



(52) 



We can now readily show that the Bogoliubov transformation in this case is 
given by 



/ \ 

a,i(t) 



y at(t) ; 



( 



2 

E 



<Xij(t) Pij{t) 



\ ( ^ 

a, 



(53) 



V 3 J 



where 



a, 



;(t) = i (B^t) - B^t) 



(54) 



and 



P ij (t)=i(B jj B ij (t)-B jj B ij (t)) 



(55) 



In terms of the above Bogoliubov coefficients the conditions (49) and (51) are 
given, respectively, by 
2 



k=i 



l 



(56) 



and 



E [onk(t)0 jk (t) - (3 ik (t)a jk (t)] = 0, 



k=l 
2 



=0. 



(57) 



k=i 
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Since Sj(t) and a\(t) are solutions of the LvN equation, the exact solution for 
the time dependent density operator with thermal equilibrium initial condi- 
tions is 



Pit) 



Z f3i Z fo 



exp < — ^ PiUJi a\(t)ai(t) + 



i=i 



(58) 



where Z^, % — 1,2, are the initial (equilibrium) partition functions for each 
oscillator as given by Eq. (17). Note that the out-of-equilibrium problem re- 
sulting from the coupling of the two systems at an initial time to is reduced 
to solving the two coupled second order differential equations given by Eqs. 
(46) and (47). As in the time dependent single harmonic oscillator case we 
can compute the time dependent thermal expectation values for the number 
operators, JVj = a\a^ 



n i {t) = Ti[p(t)a\a t ] 



i 00 l 2 r i 

y y E (ni(t),Mt)\ ex P { -HP^j ^W^'W + o 

Z /3i^/3 2 ni ,n 2 =0 { j=l 1 A 

a\a i \n 1 (t),n 2 (t)). 



(59) 



In order to evaluate the above expression we have to solve Eq. (53) for and 
al in terms of cii(t), aJt). This can be easily done using Eqs. (56) and (57) 



then obtaining 



/ \ 

a, 



\ a '/ 



2 

E 



-Pjitt) Oij(t) 



\ ( 

aj{t) 



(60) 



Using Eq. (60) in Eq. (59) we obtain 



ni{t) = E [(\<XiM 2 + \PAt)\ 2 )nj + 



(61) 
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which generalizes Eq. (25) for the case of more than one quantum oscillator. 
In Eq. (61) Uj is the initial thermal equilibrium distribution, 



n 



3 e PjUj _ i • 



(62) 



3. 1 Solution for the Instantaneous Interaction 



As a first application of the above equations, let us consider the example of two 
coupled quantum oscillators with time independent frequencies and coupling 
\{t) suddenly turned on at time t — 0, i.e, 



X(t) = 



0, t < 
A, t > 



(63) 



In this case we have for t < that 



Bij(t<0) = 



/2ui 



V 



e ZUJ 2 z 



(64) 



and for t > we can solve Eqs. (46) and (47) in terms of normal coordinates 
Qij(t), defined in terms of B^t) by 



(65) 



where the matrix transformation rjij is time independent. Substituting the 
above equation in Eqs. (46) and (47) we obtain (for details see appendix A) 



Qik(t) + n 2 k Qik(t) = 0, 



(66) 
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provided that rjij is an orthogonal matrix with elements given by 
A 2 



Vij 



1 + 



A 



(67) 
(68) 



and flj are obtained as the solutions of 



K 2 -fi 2 )K 2 -fi 2 ) = A 2 , 



(69) 



from which we obtain 



ft 



^i 2 + ^2 + VK 2 -^ 2 ) 2 + 4A 2 
! 2 + ^ 2 - V / K 2 -c 2 ) 2 + 4A 2 



(70) 



Solving Eq. (66) and substituting it in Eq. (65), we obtain for Bij(t) the 
following expression 



(71) 



k=i 



where and are constants of integration to be determined by the initial 
conditions. Since at t — 0, B^t) and B^t) are continuous, then from Eqs. 
(64) and (71), we obtain that the coefficients and 6^, appearing in Eq. 
(71), are given by 



O-ik — 

hk = 



Vik 



Vik ( 1 _ Vi_ . 



(72) 
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Substituting Eq. (72) back in Eq. (71), we obtain the result 



B^t) = £ 



VikVjk 



k=l v 



lUJi 



cos(fM) + 77- sm.(Vt k t) 



(73) 



Using Eq. (64) for at t = in Eq. (54) we obtain for the Bogoliubov 
coefficients the expressions 



(74) 



and 



2 J 



(*) + 



(*)■ 



(75) 



Eqs. (73) - (75) together with Eqs. (67) - (70) completely solve the out-of- 
equilibrium evolution of the system in consideration. In this case the system 
never reaches a final thermal equilibrium state since the Bogoliubov coef- 
ficients obtained using the result for B^t), as given by Eq. (73), have an 
oscillatory behavior for all t. The solutions for ni(t) and n 2 (^) are shown in 
Fig. 2 for the case of uoi = uo 2 . This will be further discussed below. 



3.2 An Exactly Soluble Smooth Quench 

As an example of non-instantaneous quench and exactly soluble model, let us 
consider next the case of a smooth quench with the linear coupling given by 
the expression 

A(f) = \ [1 + tanh(f/r)] , (76) 
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and with time independent frequencies U\ = uj 2 = 00. The equations for the 
Bij(t) functions now become 

B n (t) + cu 2 B, tl (t) + ^ [1 + tanh(t/r)] B l2 (t) = (77) 



and 



B i2 (t) + co 2 B i2 (t) + ^[1 + tanh^/r)]^!^) = 0. (78) 



Eqs. (77) and (78) can be more easily solved if we make use of the following 
transformation: 

R i {t)=B a {t)-B i2 {t) 

G i (t)=B il (t) + B i2 (t), (79) 

which then lead to the following uncoupled differential equations for Ri(t) and 
Gi(t), 

Ri(t) + L 2 - ^[1 +tanh(t/r)]| J R i (t) = (80) 



and 



Gi(t) + L 2 + ^[1 + tanh(t/r)]| G l (t) = 0. (81) 



Note that even for the case of time dependent frequencies u>i(t), the Eqs. (77) 
and (78) still uncouple and can be put in the form of Eqs. (80) and (81). The 
solutions for Eqs. (80) and (81) can now be found to be given by 

Ri(t) = A ie - iwt F[- ai , -61, c*, -e 2t/T ] + B^Flh, a u c, -e 2t/r ] (82) 
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and 

Gi(t) = Qe-^Fi-ai, -b 2 , c*, -e 2t ' T ] + D t e^F[b 2l a 2 , c, -e 2t % (83) 
where 



ai = y(cu - n 2 ), 
a-i = —(u - ^i), 

&2 = y(w + ^i), 

C = 1 + 

(84) 



and Qi( 2 ) = y^ 2 ± |A|. We assume that a; > |A| in order to avoid runaway 
solutions. Inverting the transformation given by (79) we now obtain that 

B a = l - {Qe-^Fi-a,, -b 2 , c*, -e 2t ^} + D t e^F[b 2 , a 2 , c, -e 2t ^} 

+ A i e-™ t F[-a 1 ,-b 1 ,c*,-e 2t / T }+ Bie^Fih,^,^-^^}} , (85) 



and 



Bin 

12 2 



\ {C ie -^F[-a 2 , -b 2 , c*, -e 2tlT \ + D^Ffa, a 2 , c, -e 2t ' T ] 

A t e- wt F[- ai , -6i, c*, -e 2 ^] - fl^Ffr, oi, c, -e 2t ' T }} , (86) 



Using the initial conditions (at t — > — oo), 
exp(^) 

-Dll = -D22 = 



B 12 = B 21 = 0, (87) 

we obtain the constants Ai,B i: Cj and A as given by: A 1 = A 2 = C\ = C 2 = 
and B\ = —B 2 = D\ = D 2 = l/y/2u. Then, finally, we find the solutions for 
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the time dependent factors Bij(t) as given by the expressions 



B n (t) = B 22 (t) = {F[h, q 1? c, -e 2t ' T ] + F[b 2 , a 2 , c, -e 2 ^]} , 

B 12 (t) = B 21 (t) = ^-j= {F[h, 01, c, -e 2t ^] - F[fc 2 , a 2 , c, -e 2 ^]} . 



The asymptotic behavior (t ^> r) for the functions Bij(t) follows from the 
asymptotic behavior of the hypergeometric functions [16], which leads to the 
results 



£n(*->oo)= ' 



r(c)r(ai - 61) ifl2t r(c)r(6i -ai) in2t 

L r (°i) r ( c - fe i) r(6i)r(c-ai) 

r( c )r(a 2 - 62) r(c)r(& 2 - a 2 ) . nit 

^ T(a 2 )T(c-b 2 ) T(b 2 )T(c-a 2 ) 



(89) 



and 



£i 2 (* -> 00) 



v/8c^ 



r(c)r( ai - ^ r^r^-aQ ^ 

r(c)r(q 2 - 62) _ T(c)T(b 2 - q 2 ) init " 

r(q 2 )r( c -6 2 ) r(6 2 )r( c -q 2 ) 



.(90) 



Note that Eqs. (89) and (90) are similar to Eq. (71) for the case of an in- 
stantaneous quench and with oj\ = u 2 . Therefore we expect that both cases 
will lead to qualitative similar results in the asymptotic limit with both os- 
cillators never reaching a final equilibrium state. This is a consequence that 
in both cases the models considered are purely reversible systems and there- 
fore energy is transferred in a reversible way from one oscillator to the other, 
with the system never actually equilibrating or thermalizing. This situation is 
changed considerably when the system is coupled to an environment, modeled 
for instance as infinity other oscillators coupled to one given oscillator, just 
like in the Caldeira-Leggett type of models [17]. Under this situation energy 
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can now be transferred irreversibly from the system of oscillators to the envi- 
ronment leading eventually to equilibration and thermalization of the system. 
How these processes come about will be explicitly shown in the next section 
in the context of the LvN approach. 

4 A system in interaction with a bath: an oscillator interacting 
with N — > oo other oscillators 

In this section we consider a system composed by an oscillator interacting lin- 
early with TV — > oo other oscillators (here to be considered as the environment 
or bath, in the typical denomination of Caldeira-Leggett type of models [17]). 
The Hamiltonian is now described by 



In the above Xk(t) is a time dependent coupling constant (with dimension of 
the square of a frequency). In general the frequencies can also be made time 
dependent, but this case will not be explicitly considered here. 

We consider that initially the oscillator and the environment are each sepa- 
rately at thermal equilibrium states at inverse temperatures f3 and f3, respec- 
tively (we do not concern here about how this initial state can be realized 
physically). That is, we consider that initially \k(t) = 0. Then the initial 
density operator p , of the total system will be given by 



1 ur 1 
H = + y ?o + 2 E [Pi + "til + 2A fc (t)g & • 



(91) 



fc=i 



Po = fa Q ® fa, 



(92) 
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where pp and pp are the initial thermal equilibrium density operators. We 
want to investigate how the initial thermal equilibrium state for the oscillator 
evolves in time when it is coupled to the environment. More explicitly, we 
want to study the equilibration process and in particular we would like to 
know how is the final equilibrium distribution reached and how long it takes 
to reach that state. 

In order to solve for the density operator in an exact way we proceed in an 
identical way as in last section. It is clear that analogous expressions to those 
obtained in the last section will also be obtained here, with similar arguments 
from the LvN approach being carried out to this problem here. For example, 
the thermal expectation value of the number operator can be shown that can 
be put in the form, in terms of Bogoliubov coefficients, as 

N 

Mt) = E [(l<v(*)l 2 + I/MOI 2 ) ^ + I/MOI 2 ] > ( 93 ) 

where p = (0,k), k = 1,2.. and n , are the initial thermal equilibrium 
distributions for the oscillator and the environment, respectively, given by 



and 



e /3u) fc _ I ' 



(95) 



Note that in Eq. (94) we have used to instead of u . As will become clearer 
in the derivation below, u>o will represent a bare frequency, while zero point 
corrections coming from the coupling to the (infinity) environment oscillators 
leads to the definition of a physical (renormalized) frequency, that here we 
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denote it by lo. The Bogoliubov coefficients are given by equations analogous 
to Eqs. (54)-(55) 

CV(*) = % (Bl v B^(t) - Bl u B, v {t)) , (96) 
M*) = * (B w B^{t) - B vv B^{t)) , (97) 

where B^ are the non zero initial values of B^ u (t), 

= 6^^=, (98) 



where = (a), a;*,). Of course, the Bogoliubov coefficients also follow identical 
relations as given by Eqs. (56)- (57). The coefficients B^t) also obey equations 
similar to Eqs. (46)-(47), 



N 



B^{t) + ulB^it) + ^k(t)BAt) = 0, (99) 



fc=i 



^(i) + ^^ fe (0 + X k (t)B^ (t) = 0. (100) 

In general Eqs. (99)- (100) cannot be solved analytically for any arbitrary time 
dependence of X(t). There is a situation however in which these equations can 
be solved exactly, and this is the case when X(t) is suddenly turned on at some 
given time, for example at t — 0. For convenience we can then express Xk(t) 
as 



X k {t) 



0, t < 

(101) 

a;fc\/^ScJ, t > 



where g is now a constant (coupling constant) parameter introduced for con- 
venience and Au; is defined below. Also in order to simplify our calculations 
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we start with uj k given by 



io k = -j-, k = 1,2, . . . 



(102) 



and then we go to the continuum limit at the end by letting L — > oo. In Eq. 
(101) Aw = f is the interval between two neighboring frequencies. For t < 
the solution of the system of Eqs. (99)-(100) is given by Eq. (98). For t > 
the solution is found by introducing, as in Sec. Ill, normal coordinates Qpuif) 
now defined by 



N 



p=0 



(103) 



where Qp P {t) satisfy 



Qw(t) + rtiQ lip (t)=0, 



(104) 



and i]^ u is an orthogonal matrix with elements (for details see appendix A) 
given by 



Von 



N 



i + E 



(105) 



(106) 



with normal frequencies obtained as the solutions of 



o _ 2 _ ^ guj 2 k Au 

fc=l u > 



(107) 



Eq. (107) can be rewritten as 



LO. 



■o 2 -^-^ = E^n2- 



(108) 
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By taking the limit iV — > oo in Eq. (108) we see that the above equation make 
no sense, since in such case the left hand side is infinity while the right hand 
is not. This is a consequence of the (infinity) zero point energy added as we 
let N — > oo. In order to overcome this problem we introduce a renormalized 
frequency uj defined as 



UJ 



= ^JujI- gNAuj. (109) 



In order to uj be finite we see that uJq have to be infinity. Then we conclude 
that ujq is not the physical frequency of the oscillator, it is a bare (unrenormal- 
ized) frequency. The physical frequency is the finite renormalized frequency 
uj. Also the introduction of the renormalized frequency guarantees that all the 
solutions VL 2 of Eq. (108) be positive allowing us to avoid possible runaway 
solutions [18]. 

Using Eq. (102) and taking N, L -> oo in Eqs. (105)-(107) we obtain 

?7o u = hm , 110 

M A "->o -cu 2 ) 2 + ttV^ 

= 2gu k Au tt^ 



and 



2VL Auj ( 2Cj 2 \ ,„ n . 
cot {LQ) = — + — 1 - — — , 112 
v 1 irg vrfi I gAuj I K ' 



where in Eq. (112) we have dropped the index ji for simplicity and it is to be 
understood that are the roots of Eq. (112). 
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Substituting the solutions of Eq. (104) in Eq. (103) leads to the following 
expression for B^it), 



N 



(113) 



p=0 



From the continuity of B^ v {t) and of its first derivative at t = we can obtain 
a^p and b w (see appendix A). Substituting the obtained expressions for these 
coefficients in Eq. (113) we obtain the result 



(114) 



from which the Bogoliubov coefficients follow 



<V(t) = \ ^B, v (t) - -^=B, v (t) 
V ^ v &*)v 



(115) 



and 



/WO = y + t^T^W- 



(116) 



Taking /x = in Eq. (93), we obtain the time dependent thermal expectation 
value of the number operator for the oscillator q Q , 



n (t) = [(\a 00 (t)\ 2 + \(3 00 (t)\ 2 ) n + \(3 o(t)\ 2 

+ E[(l«o fc (0l 2 + l«)l 2 )^ + l«0l 2 



(117) 

Then in order to obtain n (t) we have first to compute aoo(t), aok(t), Poo(t) 
and Pko(t)- These coefficients are obtained in Appendix B. They are given by 
Eqs. (B.6), (B.7), (B.14) and (B.15). From Eqs. (B.6)-(B.7) we can see that 



a 00 (t — > oo) = (3oo(t ->• oo) = 0. 



(118) 
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On the other hand from Eqs. (B.14)-(B.15) we obtain 

aot(f _ = m* : rf 2 i; 2 2a2 ::t ] a*) 



4(0-1-^)^ + ^9^1 
ft.fr .o) = J^ (g - ;f^/^7' ^.- (1 20) 



Now substituting the above equations in Eq. (117) and going to the continuum 
limit we obtain for n (t — > oo) the following result 

2g7 uj(uj 2 + uj 2 ) 1 



n (t -> oo) = /" 



[4(^2 _ £2)2 + ^2^2] ( e /3u 



,~^2 



uj J A(uj 2 - uj 2 ) 2 + ix 2 g 2 uu 2 

Note that, as a consequence of Eq. (118), the asymptotic result for the number 
density of the system's oscillator q does not depend on its initial temperature, 
only on the (continuum) ensemble initial temperature T = 1/(3. We also note 
that while the first term on the right-hand-side in Eq. (121) is convergent, 
the second term diverges logarithmically for uj — > oo. It is easy to prove that 
this divergent term is equal to the asymptotic value of (O\al(t)a (t)\0) , where 
|0) means the initial ground state of the system, i.e, the state in which all 
the harmonic oscillators are in their ground states. Then we can interpret the 
divergent term as the particles of frequency uj produced by the modification 
of the initial ground state. Clearly this term does not have a thermal origin. 
It is a response of the ground state to the interaction. In computing the final 
equilibrium temperature we have to use only the part of fi§(t — > oo) that have 
a thermal origin and this corresponds to the first (finite) term in Eq. (121). 
Alternatively note also that the last term in Eq. (121) can be eliminated by 
making a normal ordering of the initial annihilation and creation operators. 
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The distribution given by Eq. (121) is one of equilibrium, but how is it as- 
sociated to a distribution of thermal equilibrium as obtained from a thermal 
expectation value for the system's number operator ? For this, let us start by 
computing the thermal expectation value of the number operator correspond- 
ing to the particle oscillator by assuming that it is in thermal equilibrium, at 
some temperature 1//3', with the thermal bath. We then have that at thermal 
equilibrium the density operator for the total system is 



Tr(e-^) ' 



(122) 



where H is given by Eq. (91) evaluated at t — > oo, that we here express it 
directly in its diagonal form: 



N 



^ji A* + 2 



n 



(123) 



where A^ and A^ are the creation and annihilation operators associated with 
the normal coordinates that diagonalizes the original Hamiltonian and whose 
diagonalization transformation between the coordinates and normal coor- 
dinates is given by the same matrix transformation r]^, of Eqs. (110) and 

(HI), 



N 



(124) 



Using Eq. (124) it is easy to show that the annihilation and creation operators, 
a M and a*, associates to the coordinates g M , in terms of and A^ are given 
by 



N 



u=0 



\l (jj. 




(125) 
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and 



N 



a] r = V UtL 



!/=0 



(126) 



where u; M = (a), a;*;). We can now compute the interacting system's equilibrium 
number distribution, that we denote as n' , which is given by 

_ # _TV(a o a e-^) 



(127) 



Substituting Eq. (123) and Eqs. (125)-(126) with // = in Eq. (127) we get 



n 



(128) 



Finally, by taking N — > oo, L — > oo and using Eq. (110) in Eq. (128) we 
obtain, in the limit of the continuum, that 



i 



Co J [4(fi 2 - Co 2 ) 2 + n 2 g 2 Q 2 ] (e^' n - 1) 



Cj J 



VL(VL-Cuf 



A(n 2 - Co 2 ) 2 + n 2 g 2 n 2 



(129) 



Note that the above equation is identical to Eq. (121) if we identify f3' = 
(3. Then, we have showed that the particle oscillator reachs a final thermal 
equilibrium state at a temperature Tf that is exactly equal to that of the 
thermal bath, T. By introducing dimensionless parameters, 9 = % and (3 = Cuj3, 
we can write for the finite part of Eq. (121), that we denote as n eq , the following 
expression 



x(x 2 + 1) 



[4{x 2 - l) 2 + 7T 2 g 2 x 2 } (eP* - 1) 



(130) 
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In Fig. 3 we show the numerical results for n eq , normalized by a Bose-Einstein 
distribution at temperature T and energy cu, as a function of scaled temper- 
ature and coupling constants, T/uj and g/u, respectively. The results shown 
in Fig. 3 show that Eq. (130) quickly approaches to the distribution h = 
l/[exp(/3o)) — 1] at high temperatures T/uj > 1. 



5 Nonlinearities at Gaussian level 



In this section we shall consider the effect of the nonlinearities at the Gaus- 
sian approximation. The model that we consider is described by the following 
Hamiltonian with quartic potential, 

k Jt 2 + qi f + m?. (131) 

As in Sec. II we introduce time dependent annihilation and creation operators, 
a(t) = A(t)q + B(t)p, (132) 



and 



(133) 



from which using Eq. (7) we can solve for q and p in terms of a(t) and a^(t), 
obtaining 



q = i 



B*(t)a(t)-B(t)a)(t) 



(134) 
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and 



p = — i 



A*(t)a(t)- A(t)a\t) 



(135) 



Replacing Eqs. (134) and (135) in Eq. (131), normal ordering the time depen- 
dent annihilation and creation operators and truncating the Hamiltonian up 
to quadratic terms in these operators we obtain the result 



\A(t)\ 2 + ^(t)\B(t)\ 2 + ^\B(t)A a\t)a(t) 



A 2 {t) + uo 2 {t)B 2 {i) + ^&B*(t)B 3 (t) 



[A*(t)} 2 + u 2 (t)[B*(t)} 2 + ^-B(t)[B*(t)} 3 ) ^ 



. (136) 



Replacing a(t) in the LvN equation and using for H its Gaussian approxima- 
tion, as given by Eq. (136), we obtain the following equations for A(t) and 
Bit), 



A[t) = -Bit) 



(137) 



and 

Bit) + uj 2 it)BH) + ^Sl\B(t)\ 2 B(t) = 0. (138) 

Note that now Bit) is described by a nonlinear equation, that is, this shows 
how at the Gaussian level approximation it is possible to include the nonlin- 
earities of the model. As in Sec. II, the time dependent thermal average for 
the number occupation operator is given by Eq. (25), where the Bogoliubov 
coefficients are given by the same relations of the Sec. II, the only difference 
is that now Bit) is described by the nonlinear equation (138). 
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In general Eq. (138) cannot be solved exactly. But we can extract some quali- 
tative conclusions in the case in which the frequency is time independent and 
when X(t) evolves just like as in Eq. (63). In this case we have that 

Bit < 0) = -= (139) 



and for t > 0, B(t) is described by 

B{t) + u 2 B{t) + ^\B(t)\ 2 B(t) = 0. (140) 



We have to solve the above nonlinear equation with the initial condition (at 
t — 0) given from Eq. (139), 

B =vb ^ Vf (141) 



To solve Eq. (140), it is useful to express B(t) in the polar form, 

B(t) =^(t)e^ (t) , (142) 

where ip(t) and <f>(t) are real functions. Substituting Eq. (142) in Eq. (141) 

and identifying real and imaginary parts of the resultant equation we obtain 
the coupled set of differential equations 

m + (u 2 - im 2 ) m + ^\t) = o, (143) 

2i>(t)j)(t)+4>(t)4>(t)=0. (144) 

From Eq. (144) we immediatly obtain that 

0(t)=c^ 2 (t), (145) 
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where c is a constant that can be fixed by using the initial conditions. Compar- 
ing Eqs. (139) and Eqs. (142), we conclude that at t — 0, tjj — 1/V2uj, ift — 0, 
= and = uj. Using these values in Eq. (145) we obtain that c = 1/2. 
Then integrating Eq. (145) with the initial condition = we obtain, 

1 r 

4>{t) = -Jdt'r 2 {t'). (M6) 



To solve for ip(t) we replace Eq. (145) (with c = 1/2) in Eq. (143) to obtain 
the result 

m + - \r\t) + ^ 3 (t) = o. (147) 



Note that the above equation is singular for -0 = 0, but as we shall show below, 
this singularity is never reached by ip(t) in the course of its time evolution. 
Using Eqs. (21) and (22) in Eq. (25) and replacing B(t) as given by Eq. (142), 
we obtain for the time dependent thermal expectation value of the number 
operator, n(t), the following expression 



fi(t) = 



i> 2 (t)+u; 2 i> 2 (t) + -r 2 (t) 



n 

UJ 



+ 



1 



1 

2" 



(148) 



Then, in order to know the time evolution of n(t) we only need to know the 
time evolution of ip(t). From Eq. (147) and using the initial conditions for ip(t) 
and ip{t) we obtain 



(149) 
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from which we get 



dx 



t = ± 



/ 




A 



) 



X — 0J 2 X 2 — jX 3 



(150) 



l/2w 



16w 2 



The above integral cannot be solved analytically, and then we are not able 
to solve for ip{t). But we can obtain some qualitative conclusions about the 
behavior in time of ip(t) that will help us to understand the effect of nonlin- 
earities on the thermal equilibration of the system we are studying. For this 
end we rewrite Eq. (148) by using Eq. (149), obtaining 



Then, if we know the qualitative behavior of ip(t) we would equally well know 
the qualitative behavior of n(t). If in Eq. (149) we think in ip{t) as being the 
analogous of the coordinate of a particle, then Eq. (149) is nothing but the 
energy conservation law of a particle of unity mass that moves in a conservative 
potential U (ijj) given by 



It is clear from the potential U (ip) that for any finite energy ip(t) evolve in time 
in an oscillatory way, never reaching a fixed value, and therefore the system 
never reaches a final thermal equilibrium state. Note also that in order to ip(t) 
to reach the singularity — 0, the energy, in the particle mechanical analogy, 
need to be infinity. Since the energy, as given by the right hand side of Eq. 
(149), is finite, we also conclude that such singularity can never be reached. 




(151) 



(152) 
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6 Conclusions 



In this work we have studied several simple nonequilibrium models represented 
by quantum oscillators in interaction that can be exactly solved. In special we 
have studied examples of smooth and sudden quenches that were explicitly 
solved in the context of the LvN approach. As emphasized in Refs. [11,12], 
the LvN approach has a clear advantage over other techniques to study out-of- 
thermal equilibrium problems since standard methods of quantum mechanics 
are used. We have explicitly seen this here by working directly in terms of 
annihilation and creation operators and their Bogoliubov transformation re- 
lations, allowing to determine the evolution of number densities among other 
quantities, like the final equilibration temperatures and the general system 
time evolution. 

Previous use of the method, like to the study of second order phase transi- 
tions [12], for the decoherence problem [19] and for the studies performed here, 
show that the LvN approach can easily be adapted to study more complex 
problems that are difficult to deal with in the context of nonequilibrium dy- 
namics of fields where usually the expressions for nonequilibrium propagators 
and correlation functions must be worked out. In this case the LvN approach 
may provide an avenue to more rigorously take into account in the dynamics 
of systems the interactions with the environment degrees of freedom. An ex- 
ample of this has been worked out in details in Sec. IV, for the model of the 
system quantum oscillator in interaction with an infinity set of other quantum 
oscillators. 

The extension of the method for applications to the current important prob- 
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lems of decoherence and description of dissipative dynamics of quantum sys- 
tems, for the more realistic cases where nonlinear interactions are present, is 
currently in progress and we expect to report on them in the near future. 
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A 

In this appendix we solve the system of equations given by 

N 

B^{t) + tv 2 B^(t) + hB^t) = 0, (A.l) 



When N — 1, k — 1 and then we can write A^ = A the equations above 
become identical to Eqs. (46)-(47) with time independent coupling constant 
(and changing the indexes 0, 1 to 1, 2). On the other hand, when A& = Uk^gAu 
we recover Eqs. (99)-(100) with time independent coupling constant as given 
by Eq. (101). Also in this last case, in the initial conditions given by Eq. 
(A. 3) we have to replace the bare frequencies by the renormalized ones 



B^it) + ulB^it) + \ k B^(t) = 0, 



(A.2) 



with initial conditions at t — given by 




(A.3) 
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oj^ = (iv,Uk) [see Eq. (98)]. In order to solve Eqs. (A.1)-(A.2) we introduce 
normal coordinates Qpvit) through 



N 



p=0 



(A.4) 



where only the second index of B^ u (t) appear in the matrix transformation 
{rj} since the first index are uncoupled in Eqs. (A.1)-(A.2). Substituting Eq. 
(A.4) in Eqs. (A.1)-(A.2) we obtain 



N 



k=l 



N 

£ 

p=0 
N 

£ [vkpQpp(t) + ^hk P Qpp{t) + XkVopQpp(t) 

p=0 



o, 



= 0. 



(A.5) 



(A.6) 



Multiplying Eq. (A.5) by r] 0cr and Eq. (A.6) by r] ka and adding the resulting 
equations we obtain 



N 

£ 

p=0 



N 



N 



J2 t^lVvpVva + J2^k (VkpVOa + VOpVka] 



.u=0 



k=l 



N 



+ J2v^v^Qpp(t) \=o, 



and imposing that, 



N 



^ ' TjvpT]ua $pcri 



u=0 



N N 

£ ^IVupVua + £ Afe (l] kp 1] 0r7 + 1] 0p 1] ka ) = S p(T fll, 
u=0 k=l 



(A.7) 



(A.8) 
(A.9) 



we obtain 



Qp,(t) + nlQ^(t) = o. 



(A.10) 
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Eq. (A. 9) can also be written as 

TV 

E 

fc=l 



N 

E 

u=0 



1]ua ^per^cr- 



(A.11) 



From Eqs. (A. 8) and (A. 11) we obtain that 



N 

^IVup + E ^ (Vkp$ov + Vophu) = V»ptt 
k=\ 



(A.12) 



Taking u = in the above equation we obtain that 



N 



CO, 



oVop + E X kVkp = VOp^p, 

k=l 



(A.13) 



while for v ^ we obtain for r\ kp the result 



Vkp = 



(A.14) 



Taking p = a in Eq. (A. 8) and using Eq. (A.14) we can solve for i] 0p obtaining 
the expression 



N \2 

i + e 



(A.15) 



Using Eq. (A.14) in Eq. (A.13) we also get that 

N 



N \2 
Ai, 



k=l U k 



(A.16) 



Now, from Eq. (A. 10) we can write for B^ v {t) the following expression 



N 



p=0 



(A.17) 



and using the initial conditions given by Eq. (A. 3) we get for a w and b pp , 
respectively, the equations 
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v n„ 





(A.18) 



Substituting Eq. (A.18) in Eq. (A. 17) we finally obtain the result 



N 



B,At) = E ^ 



CO 



l + zJL) e ^pt + I i _ _/± ] e 



CO 



a 



(A.19) 



In this section we deduce the Bogoliubov coefficients aoo(t), Pooit), aotit), and 
Pko(t). Using Eq. (114) in Eqs. (115)-(116) we obtain 



N 



a 



(*)=E 



p=0 k V 4fi P 



(Q p + u^) e 1 ^ + (^ - c^) e^ 



, (B.1) 



TV 



p=0 k v 4fi P 

VupVvp 



u a + a,) e if ^ + (u) u - n o ) e 



Taking (j,, v — in the above equations we obtain 



"oo = 2^ —r 

p=0 * 



2 + — + ^ e inpt + I 2 - — - ^ I e 



-iQ„t 



Q p Co 



Q p co 



A» = E 



V (Vop) 



p=0 



Qp to J \uj Q P J 



(B.2) 



(B.3) 
(B.4) 



Substituting Eq. (110) in the above equations and going to the continuum 
limit we obtain, e.g. for aoo(t), the expression 
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«ooW = % I dVL- 
00 J 





n 



2„-iClt 



4 (Q 2 - Co 2 )' + 7i 2 g 2 Q 2 



JL 

2Co 
9_ 

00 



dVL- 



(Co + ft)V m -(Co- Vl) 2 e 



— oo 
oc 



4(n 2 - co 2 ) 2 + n 2 g 2 n 2 
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The integral in Eq. (B.5) can be easily evaluated using the Cauchy theorem and 
choosing as contour of integration in the upper half plane. The poles located 
in the upper half plane are ±u + where c<j = Jui 2 — ZE jf-. Assuming a; > ^ 
the result is 
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In the same way we obtain from Eq. (B.4) the expression for fioo(t) (assuming 



again u > 



nge 



-■ngt/A 



ixg - 4iu?j e llJjt - (ng + 4iio^ e 



—iuit 



(B.7) 



Setting /i — 0, v — k in Eq. (B.l) we obtain, 
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4 (Q 2 - a; 2 )^ + ttV^ 2 
In going to the continuum limit we have to take care since there will appear 
a singularity when Q = io k . In reality this singularity is fictitious and it can 
be eliminated from Eq. (B.8) using Eq. (108). In terms of the renormalized 
frequency Co Eq. (108) can be written as (setting ji — k) 
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where we have isolated the pole at Uk = &k and the prime in the integral means 
that the point uj = Qk has been excluded from the integration. It is easy to 
show that such integration is zero. For this end we write such integration as 



(B.10) 



where now it is understood that we omit the points ±Qk = ^k in doing the 
integration over the frequency. Going to the complex plane and using Cauchy 
theorem we obtain, 
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Then, from Eq. (B.9) we obtain 
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where it is understood that Q k = ^k and since we are considering the contin- 
uum limit, in which case Alu — > 0, then the right-hand-side of Eq. (B.12) is 
well defined. Now we can use Eq. (B.12) in Eq. (B.8) to eliminate the singu- 
larity after isolating it. Isolating the singularity in Eq. (B.8) and going to the 
continuum limit we obtain 
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where in passing to the second equality we have used the even symmetry of 
the integrand, factored the term (oo k + fi), also we have used Eq. (B.12) to 
eliminate the singularity (Q k — oo 2 )~ l and finally we set Q k = oo k . The prime 
in the integrand of Eq. (B.13) means that we are avoiding the points ±uo k . 
Then when using the Cauchy theorem we have to multiply the residues in 
these poles by in instead of 2m. Actually the pole at — uj k has been factored. 
Then after using Cauchy theorem it is easy to obtain, 
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Setting \i = k, v = in Eq. (B.2) we obtain an expression for /3 k o(t) that 
contains also the same singularity as in the case of a ok (t). Then, proceeding 
in the same way as in the above case we obtain 
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B.l. The number density for ujq/3 = 1 and for decay times (in units of 1/u>q) 
1.0, 1.5 and 2.0. 
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Fig. B.3. The system's equilibrium distribution, Eq. (130), normalized by the dis- 
tribution no = l/[exp(/3u>) — 1], as a function of the (scaled) temperature T and for 
different values of the coupling constant g. 
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